Phantom Field with O(N) symmetry in Exponential Potential 
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In this paper, we study the phase space of phantom model with O(N) symmetry in exponential 
potential. Different from the model without O(N) symmetry, the introduction of the symmetry leads 
to a lower bound w > —3 on the equation of state for the existence of stable phantom dominated 
attractor phase. The reconstruction relation between the potential of O(N) phantom system and 
red shift has been derived. 
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I. INTRODUCTION 

Astronomical measurements from Supernovae [2- 13 
and CMB anisotropy 0,0,0 independently confirm that 
about two-thirds of the energy density in our Universe is 
dark energy whose effective equation of state may lies in 
the range —1.38 < w < —0.820- A recent analysis for 
the new high red shift dataQ as well as dark energy and 
dark matter halos[? ] also indicate that dark energy with 
w < — 1 seems more favorable. While the equation of 
state of conventional quintessence models lOjthat based 
on a scalar field and positive kinetic energy can not evolve 
i<> the tin- regime of <r — 1. Some mil hors |y||^|^|^ 
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models which have negative kinetic energy and can realize 
the w < —1 in its evolution. Although the introduction 
of a phantom field causes many theoretical problems such 
as the violation of some widely accepted energy condition 
and the rapid vacuum decay |2lJ, it is still very interesting 
in term of fitting current observations. Comparing with 
the other approaches to realizing the w < — 1, such as the 
modification of Friedman equation, it seems economical. 

On the other hand, the role of complex scalar field as 
quintessence has been studied in0 an d its generaliza- 
tion to the fields with O(N) symmetry was done in Ref. 
|4l| . In this paper, we study the case that when phantom 
models have an O(N) symmetry constraint via a specific 
model in which the potential is exponentially dependent 
on the field. We show that the dynamical system admits 
only the phantom dominated attractor phase and the in- 
troduction of O(N) symmetry leads to a lower bound 
on the equation of state as w > —3, which is absent 
in the singlet scalar field phantom models. It is worth 
noting that although the O(N) phantom and the O(N) 
quintessence lead to similar equations, but they have very 
different dynamical evolution and physical implications, 
especially they will lead to different evolution of the equa- 
tion of state w. We also derive the relation between the 
red-shift and the potential of the scalar fields, known as 



reconstruction relation. It turns out that the potential 
relates to the red-shift the same way as that in ordinary 
scalar fields theory while the variation of the fields does 
in a quite different manner. 



II. O(N) PHANTOM 

We start from the flat Robertson- Walker metric 



ds 2 = dt 2 -a 2 {t)dx 2 



(1) 



The Lagrangian density for the Phantom with O(N) sym- 
metry is 

L 9 = -\g^(d^ a ){d^ a ) - U(|$ a |) (2) 

where $ a is the component of the scalar field, a = 
1,2,- • • , N. To make it possess a O(N) symmetry, we 
write it in the following form 



$ x = R(t) cospi(f) 

$ 2 = R(i) sin tpi (t) cos <p 2 (t) 

$ 3 = R(t) sin ipi (t) sin <p 2 (t) cos (fs (t) 



(3) 



N-l 



= R(t) sin <px (t) • • • sin< / 97v-2(i) cos (fN-x(t) 
<& N — R(t) smipx(t) ■ ■ ■ sin^ A r_ 2 (£) sinipiy^i(t) 



Therefore, we have |$ a | — R and assume that the poten- 
tial of the O(N) phantom depends only on R. 
The action for the universe is : 



S 
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(4) 



'Electronic address: kychz@shtu.edu.cn 



where g is the determinant of the metric tensor g^ LV , G 
is the Newton's constant, R s is the Ricci scalar, and p 7 
is the pressure of the baryotropic fluid whose equation 
of state is p 7 = (7 — l)p 7 . Varying the action, one can 
obtain the Einstein equations and the radial equation of 
scalar fields as 
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V(R)} 



with p$ = -\{R 2 + ^) + V(R) and = -±(i? 2 + 

Jfjp ) — V^-R) are the energy density and pressure of the 
$ field respectively, and H is Hubble parameter. To ob- 
tain the above equations, we make use of the fact that 
the " angular components" of the system could be simpli- 
fied by the first integrals and its net contribution to the 
dynamics could be manifested by the term containing f2 
in the radial equation. For details of this, one can refer 
to Ref . [4l| . The equation of state for the O(N) phantom 



dynamical system. The variables could be defined as 

TV = logo. The equation svstem(ISll- Hllf) become the fol- 
lowing autonomous system: 



dx 
dN 



dy 
dN 
dz 

In 

dN 



-£[7(1 + x 2 - y 2 + z 2 ) - 2{x 2 + z 2 ^ 



(12) 



-(3 x -V6z 2 t+ x /-\y 2 ) 



'-y[j(l + x 2 - y 2 + z 2 ) - 2(x 2 + z 2 ^ 



Xxy 



-3z + -z[j(l + x 2 -y 2 + z 2 ) - 2(x 2 + z 2 )] - V6xz£, 



Also, we have a constraint equation 



w 
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It is clear that the O(N) phantom could realize the 
equation of state w < — 1 which is equivalent to 



< R z 



n 2 



< 2V(R) 



(7) 
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where the term \ Rl comes from the "total angular mo- 
tion". 



III. ATTRACTOR PROPERTY OF O(N) 
PHANTOM 



In this section, we firstly investigate the attractor prop- 
erty of the O(N) phantom field in an exponential poten- 
tial. To do so, we rewrite the equations of motion as 



H 



- (pj +p~ l - R 2 
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where 



R 3H 2 



n K PR 2 2 2 



3H 2 



(13) 



(14) 



The equation of state for the O(JV) phantom could be 
expressed in term of the new variables as 



1H> 
P<1> 



x 2 — y 2 



(15) 



The critical points of the above autonomous system 
are easily obtained by setting the right hand sides of the 
above equations to zero. We write the variables near 
the critical points (x c ,y c , z c ,£c) in the form x = x c + u, 
y = y c + v, z — z c + w and £ = £ c + x, where u, v, W, x 
are perturbations of the variables near the critical points 
and form a column vector denoted as U . Substitute the 
above expression into the autonomous system (|12fl . one 
can obtain the equations for the perturbations up to the 
first order as: 



p 7 = -3H(p 1 +p 7 ) 



(9) 



U' = M-U 



(16) 
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(10) 



) + V(R)] (11) 



The potential V(R) here is specified as V(R) = 
Vq exp(— XnR). Now, we will introduce the following vari- 
ables to obtain the autonomous system for the above 



where the prime denotes the derivative with respect to N. 
The coefficients of the perturbation equations form a 4 x 4 
matrix M whose eigenvalues determine the type and sta- 
bility of the critical points. The only physically meaning- 
ful critical point corresponding to the autonomous sys- 
tem (H2J) is (x,y, z,£) = (-^7=, + ^,0,0), which cor- 
responds to the eigenvalues (0,-3 — 4r, — 37 — A 2 ,— 3 + 
4-). Therefore, it is a stable node of the autonomous 
system when A 2 < 6. This corresponds to a late time 
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FIG. 1: The phase diagram of the 0(A) Phantom system in 
terms of x, y, z for different initial x, y, z and f. 
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FIG. 2: The projection of the phase space in x y plane. 



attractor solution which is a phantom dominated epoch 
£Ir = 1 and an equation of state w = — ^ — 1. Unlike 



the singlet phantom field in exponential potential [23j, 
the introduction of internal symmetry imposes a lower 
bound to the parameter A 2 for attractor solution, which 
equivalent to the equation of state w > —3. This is a very 
interesting characteristic of the O(A) phantom field. In 
Fig. 1 and Fig. 2, we show the numerical results of the 
phase space of the system. In this plots, we choose the 
parameter 7 = 1 and A = 1.2. 



IV. RECONSTRUCTION AND DISCUSSION 

Now, we will correlate the potential with red shift. To 
do so, following the earlier study in this field^^i we m_ 
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FIG. 3: The evolution of the equation of state w vs. N for 
different initial x, y, z and f. 



troduce the quantity 



r(z) = 



to 



t(z) a(u) 7 H(x) 



dx 



(17) 



which is the Robertson- Walker coordinate distance to an 
object at red-shift z. Also, we denote 



Pm 



3r> M ffg(l + z) 3 
8ttG 



(18) 



where H is the present Hubble constant, Om is the frac- 
tion of non-relativistic matter to the critical density. We 
then readily have 



H{zf 
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(dr/dz) 2 
+ (l + z) 



(19) 



d 2 r/dz 2 

a (dr/dz) 2 ' v ' '(dr/dz) 3 

§ = -(l + z) H(z) = -(l + z) d f 
dt dz 

From the above equations, it is not difficult to derive the 
reconstruction equations as 



V[R(z)] = 



1 
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(21) 

(l + z)(d 2 /dz 2 ) 3 3 
(dr/dz) 3 + 2 [ + ' 



Eq.(Hl) is the same as those ordinary quintessence while 
the Ea. (|21(l is different in that there is a sign difference: 
the right hand side of the Ea. (|21|l is positive while it is 
negative in conventional O(A) quintessence model. 
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Up to now, the observation data do not tell us what 
should be the nature of dark energy. But the future ob- 
servation will be helpful to determine whether the dark 
energy is phantom, quintessence, or cosmological con- 
stant. If the equation of state w < —1 is completely 
confirmed by observations, then its implications to fun- 
damental physics would be astounding, since it cannot 
be achieved with substance with canonical lagrangian. 
Phantom field could be a good candidate for such sub- 
stance because of its negative kinetic energy and sim- 
plicity. In this paper, we study a class of new phan- 
tom models, in which the scalar field possesses a O(N) 
internal symmetry. In a specific model, in which the 
potential of the phantom fields is an exponential po- 
tential, we show that there exists an attractor solution 
when A 2 < 6 and the attractor solution corresponds to 
a phantom energy dominant phase and an equation of 



state w — — 4j — 1. Most strikingly, the introduction of 
O(N) symmetry imposes a restriction to the existence of 
the attractor solution A 2 < 6, which accordingly puts a 
lower bound to the equation of state w > — 3. On the 
other hand, the speed of sound of the dark energy is de- 
fined by c 2 = &i = , , where p = LM, X) 

&ndp = 2XL^ x -L 4> {(t> 1 X) with X = ^(c^) 2 . So, The 
model in this paper will produce an acceptable sound of 
speed c 2 = 1. 
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